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Abstract. We recently suggested that the Scissors Mode (a collective excitation in which one system rotates 
with respect another one conserving its shape) can occur in crystals with axially symmetric atoms as a 
precession of these atoms around the anisotropy axis of their cells, giving rise to a form of dichroism. In 
the present paper we investigate how the Scissors Mode can be realized in crystals with cubic symmetry 
and evaluate its photo-absorption cross-section. This turns out to be of the same order of magnitude as 
that for crystals with axially symmetric atoms, but does not exhibit any correlation between the direction 
of the photon and the axes of the cell. 
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1 Introduction 

The Scissors Mode is a collective excitation in which two 
systems move with respect to each other conserving their 
shape. It was first predicted to occur in deformed atomic 
nuclei [1] by a semiclassical calculation in which protons 
and neutrons were assumed to form two distinct bodies 
to be identified with scissors blades. Their relative rota- 
tion generates a magnetic dipole moment whose coupling 
with the electromagnetic field provides the signature of 
the mode. 

After its discovery [2] in a rare earth nucleus, ^^^Gd, 
and its systematic experimental and theoretical investi- 
gation [S] in all deformed atomic nuclei, it was predicted 
to occur in several other systems including metal clus- 
ters [1], quantum dots [S], Bose-Einstein [6j and Fermi [7] 
condensates and crystals |3- In all these other systems 
one of the blades of the scissors must be identified with a 
moving cloud of particles (electrons in metal clusters and 
quantum dots, atoms in Bose-Einstein and Fermi conden- 
sates, an atom in a cell in crystals) and the other one with 
a structure at rest (the trap in Bose-Einstein and Fermi 
condensates, the lattice in metal clusters, quantum dots 
and crystals). Among these other systems until now the 
existence of the Scissors mode has been confirmed only in 
Bose-Einstein condensates |9j. 

The case of crystals has been studied only for atoms 
with axial symmetry and ionic bond, in a semiclassical 
model [8| in which an atom whose charge density profile 
is approximately ellipsoidal is regarded as a rigid body 
whose symmetry axis can precess around the axes of its 
cell under the electrostatic force generated by the ligands. 
If the deformation is sufficiently large it can rotate only 
around the anisotropy axis of the cell because the restor- 



ing force is too strong to allow a precession around the 
other axes. With such a geometry the Scissors Mode gives 
rise to a magnetic dichroism, with a strong Ml transition 
amplitude which is maximum when the momentum of the 
absorbed photon is parallel to the anisotropy axis of the 
cell and zero when it is orthogonal to it. A numerical esti- 
mate was done for LaMnOs. It was found that the energy 
of the Scissors Mode is about 4 eV (but see Section 3) 
and the Ml transition amplitude of the order of 1.5 a.u.. 
For small deformations one should expect, in analogy with 
atomic nuclei [10], a splitting of this collective mode cor- 
responding to precessions around all the cell axes. 

There is a great deal of compounds whose structure 
has cubic symmetry, either at a global level (cubic space 
groups) or at a local level. We have in mind the transition 
metal monoxides that crystallize in the rocksalt structure, 
like TiO and NiO, other that adopt the rare stechiometry 
MO3, like ReOa, all the ternary compound having a per- 
ovskite cubic structure, named after the mineral CaTiOa, 
etc., not to mention the many transition metal complexes 
in octrahedral coordination, like Fe(CN)g~ or Co(NH3)g~'~. 
Among the Rare Earth compounds we mention Europium 
oxide (EuO) as a prototypical example of cubic crystal, 
but one can find other examples in cubic Laves phases of 
Rare Earth and pnictides. Examples of Rare Earth ions 
in octahedral coordination can be found also in [TT] . 

Obviously the analysis of Ref. [5] cannot be valid for 
such systems, but the magnetic nature of the Scissors 
Mode makes its investigation more interesting because 
magnetic anisotropy is at the origin of interesting techno- 
logical applications of some such systems, including mag- 
netic storage devices and sensors, spin-torque nano-oscillators 
for high-speed spintronics and spin-optics [H]. We would 
like in fact to remind that in the cases studied so far, 
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the Scissors Modes provides specific pieces of information 
about the structure of the systems. In nuclear physics it 
is related to the superfluidity of deformed nuclei, in Bose- 
Einstein Condensates it provides a signature of superflu- 
idity, in metal clusters it is predicted to be responsible 
for paramagnetism and in crystals with axially symmetric 
atoms of a form of magnetic dichroism. In a separate work 
we will relate Scissors Modes in crystals with cubic sym- 
metry in the Rare Earth region to Spin-Orbit Locking, 
a conjectured atomic structure in which, due to a strong 
spin-orbit force, the charge profile of an atom in a cell is 
firmly bound to the atomic spin. More generally the study 
of the possible occurrence of such collective modes is in- 
teresting in understanding to which extent it is a universal 
feature of finite many-body systems. 

In the present paper we investigate in which way the 
Scissors mode can be realized in crystals with cubic sym- 
metry. We find that if it is excited by a photon whose 
momentum is parallel to one of the symmetry axes, the 
motion resembles that of a balance-wheel but, unlike the 
case of atoms with axial symmetry, the absorption cross- 
section does not depend on the direction of the photon 
momentum. 

In Section 2 we describe the semiclassical model which 
we adopt and we evaluate excitation energy and photoab- 
sorption cross-section. In Section 3 we discuss the numer- 
ical estimates of the parameters of the model and their 
relation to some conceptual issues. We conclude in Sec- 
tion 4. 



2 The model 

We assume a frame of reference with the Xi, 0:2, Xa-axes 
parallel to the axes of the cell, and we denote by ^1,^2,^3 
the axes of the principal frame of inertia of the atom, 
regarded as a rigid body of cubic profile, as shown in Fig. 
1. If we call (f>i the rotation angle around the the 
Hamiltonian of the atom is 



H 
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21 



(1) 



where I is the moment of inertia of the atom with respect 
to any of its axes. This expression is justified for the lowest 
lying excitations for which, due to the cubic symmetry, 
we do not need to express the kinetic energy in terms of 
the components of the full angular momentum. Moreover, 
because of this symmetry the potential must satisfy the 
relation Ri{^)V{(j)i) — V{(f>i), where Ri{4>i) is the rotation 
operator through the angle 0^ around the i-axis. 

Therefore it is sufRcient to define it in the interval 
~ f 5[ < f , hi which we assume the harmonic form 

V{cj^) = lccb\ -\<cp<-^. (2) 

The eigenf unctions normalized in this interval for each of 
the angles (pi are harmonic oscillator wave functions 



A„(0) 



1 



v/\/^2"n! (/)( 



exp 



(3) 




Fig. 1. A section of an atom with cubic symmetry in the 
plane xi,X2- Its principal frame of inertia, with axes ^1,^2, 
oscillates around the x^, = ^3 axis as a balance-wheel. Its profile 
can approximately be obtained by a superposition of ellipsoids. 
Only the charge external to the sphere contributes to both the 
restoring force constant and the moment of inertia. 



where _ff„ are Hermite polynomials and 
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with eigenvalues 
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Because of the cubic symmetry the eigenf unctions in the 
full range of </> are 
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We will consider only transitions to the lowest excited 
states, denoted by — 1,2,3 if the oscillation takes 
place around the i-th axis, and we well denote by tfo the 
ground state. For instance 

!Z'l('/>l,'/>2>3) = V'l(0l)V'o('/'2)V'o('/>3)- (7) 

The interaction Hamiltonian with the electromagnetic field 



Hi = V --p« • A(r«) + ^ S« • V« x A(r«) (8) 



where to, e,/i, S,p and r are the mass, electric charge, 
magnetic moment, spin, momentum and position vector 
of the electron, and the sum extends over the electrons of 
the atom. Writing the vector potential as 



A ^ e'^) exp(ik • r) 



(9) 
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where e*-^) is the photon polarization vector, the absorp- 
tion cross-section for a photon of polarization A and energy 
h(jj = hkc, c being the velocity of light is 



(A) ,2 



(10) 



where a is the fine structure constant and 
Aff ^ = {W,\e'^-''e'-^^-J2 (vC) -iS(') x V^'^) \Wo) . (11) 

In our semiclassical approximation the deformed atoms 
rotate as rigid bodies. Therefore when they perform rota- 
tional oscillations around the i-th axis the velocity of an 
electron at position r, is 



dt 



(12) 



where is the unit vector in the i-th direction. The time 
derivative of (j) is determined by the equation of motion 



dt 



so that 



(13) 
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In the electric dipole approximation (neglecting the expo- 
nential ) the matrix element (111 vanishes, because the 
contributions from the four regions of (pi cancel out with 
each other. Enforcing the cubic symmetry is essential to 
get this result, which corresponds to the physical fact that 
the electron velocities are opposite in opposite regions of 
(j). Taking into account the first term in the expansion of 
the exponential we get 
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Now we make the approximation 



^m((7 



-I. 



(15) 



(16) 



Concerning the sum over the spins we observe that in 
some of the compounds with cubic symmetry the spins of 
localized electrons are not balanced. The atoms of these 
compounds have a net spin of the order of several units 
of electronic spin. If the spin-orbit force is strong with 
respect to the crystalline field, the atomic spin is firmly 
coupled to the charge density profile [13] , a coupling which 
we call Spin-orbit Locking. In this case the atomic spin 
will follow the rotation of the charge profile. Since the 
rotation angle is very small (see below), we can make the 
approximation 

^S(')«i,S (17) 
I 



where ry is the ratio of the number of electrons which con- 
tribute to the spin over the number of electrons which 
participate to the collective motion and S is the unit vec- 
tor in the direction of the spin. At last 
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We see that the contribution of the spin is negligible, so 
that study of excitation of the Scissors Modes by pho- 
toabsorption cannot give information about the Spin-orbit 
Locking structure. Another way of investigating such struc- 
ture by exciting Scissors Modes will be presented in a 
separate paper [T4j. Using Eq.Q we find the absorption 
cross-section 

--"J^^^VXC. (20) 

We can repeat the above calculation with the wave 
functions appropriate to atoms with axial symmetry [S] . In 
this case the Scissors states are labeled by the component 
of the angular momentum along the cell anisotropy axis, 
m = ±1, and we find 
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(21) 



where k, a are the unit vectors in the direction of the 
photon momentum and the anisotropy axis respectively. 
The cross-section for photon absorption by crystals with 
axially symmetric atoms is 



(22) 



It gets its maximum when the photon momentum is par- 
allel to the anisotropy axis of the cell, and vanishes when 
it is perpendicular to it. This dichroism is due to the fact 
that the angular momentum of the photon has the di- 
rection of its momentum. Since the atom can rotate only 
around the anisotropy axis of the cell, it can absorb the 
photon angular momentum only in this direction. 



3 Moment of inertia and restoring force 
constant 

In this Section we make a few considerations about how to 
make the estimate of moment of inertia and restoring force 
constant a bit more realistic remaining in the framework 
of a semiclassical model. 

In the case of atoms with axial symmetry and ionic 
bond an analytic evaluation of these parameters was pos- 
sible in the framework of the present semiclassical model 
under the assumption that the precessing atom can be 
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treated as a rigid body with constant charge density. Some 
improvement can easily be obtained by replacing the con- 
stant density by some theoretical or phenomenological 
function p{r). The expressions of moment of inertia and 
restoring force constant can no longer be given in closed 
form, but can easily be evaluated. For instance the poten- 
tial energy in the crystalline field Vc of a charge density p 
rotated through the angle around the z-axis is given by 
the expression J (irV^(r)e~*'^^'/Ci(r), where Li is the i-th 
component of the angular momentum, so that 



inertia the square of the Scissors Mode frequency becomes 



J dvV,{v)Ll p(v) . 



(23) 



More important changes can be devised thinking about 
a microscopic calculation based on a Random Phase Ap- 
proximation. In such a framework one can neglect elec- 
trons of the atomic core weekly coupled to the electrons 
responsible for deformation, which are the ones partici- 
pating to the collective motion. This feature is already 
present in our evaluation of the restoring force constant (in 
which only the electrons external to the maximal sphere 
inscribed in the charge profile contribute), but not in the 
evaluation of the moment of inertia. It can be interesting 
to compare with the corresponding situation in Nuclear 
Physics. For deformed atomic nuclei it is known that the 
actual value of the moment of inertia is smaller than that 
of a rigid body. A semiquantitative estimate in this case 
can be obtained by neglecting the contribution of the nu- 
cleons internal to the largest sphere inscribed in the nu- 
cleus as shown in Fig. 2, getting 



effective 



• M 



M 
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where X^^Yas,s\c classical value of the moment of iner- 

tia of an ellipsoid whose major and smaller axes, rMj^rm 
are parametrized according to vm — ?'o(l + 2(5), r,„ = 
(1 — 6). This has a suggestive interpretation in the fact 
that since deformed atomic nuclei are superfluid, their ex- 
ternal part can rotate with respect to the internal spheri- 
cal part without friction [TS]. But the actual justification 
is provided by the agreement with microscopic calcula- 
tions [5] . In the case of deformed atoms we cannot invoke 
the same argument, but we can ignore the electrons of the 
core which are weekly coupled to the electrons responsible 
for deformation. What is the fraction of these electrons is a 
dynamical problem. Probably the actual value of the mo- 
ment of inertia lies between the classical and the effective 
ones. 

The procedure just outlined might be regarded as a 
way to implement at the semiclassical level the quantum 
mechanical requirement that a spherical object cannot 
have a rotational energy: if the moment of inertia accord- 
ing to Eq. (24) is proportional to \6\, the rotational energy 
of the ellipsoid becomes infinitely large for \S\ — ^ and 
therefore rotational states do not exist in this limit. 

Using the expression of the restoring force constant 
determined in [8] , and the effective value of the moment of 
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where Zi is the number of ligands and Rq the distance 
of the ligands from the center of the cell. We see that 
the excitation energy does not depend on the size or the 
number of electrons of the atom at the center of the cell, 
but only on its deformation and the number of ligands, 
their charge and distance from the cell center. The re- 
placement of 2(,lassic -^effective increases the energy of 
the Scissors Mode which in the case of LaMnOs becomes 
approximately 9 eV. 

The above results can immediately be extended to 
some of the crystals with cubic symmetry. They are the 
crystals whose atoms have a charge profile which can be 
approximated by superimposing three ellipsoids with the 
same center and the axes in the directions (see for in- 
stance Ref. [1^). Therefore both moment of inertia and 
restoring force constants have twice the values for an el- 
lipsoidal atom. 

We distinguish two cases. The first one is epitomized 
by transition metal oxides of cubic symmetry (like NiO) , 
where the valence shell contributing to the mode is the 
outermost shell (3d orbitals in the first transition metal 
series). In this case the Crystal Field is of the order of 
1^2 eV, indicating that a hybridization model approach 
is very appropriate. The second one is exemplified by rare 
earth oxides of the same symmetry (like EuO), where the 
active shell (4f shell) is shielded by occupied external shells 
(e.g 5s and 5p). In this latter case the effective Crystal 
Field is very small (roughly one tenth of that of the tran- 
sition metals), which also indicates that a Crystal Field 
model might be appropriate. However this does not mean 
that the 4f electrons are free to rotate (i.e. they are not 
core states) since there is experimental evidence [12] that 
these orbitals are somehow bound to the lattice, though 
by an hybridization energy which is one order of magni- 
tude less than in the case of transition metals. Therefore 
in both cases, in order to calculate the restoring force, 
we can adopt a hybridization model. In this model, it is 
known that the hybridization energy of a bond is given 
by the gain in electrostatic energy of the electron density 
accumulated along the bond, in order to take advantage 
of the attraction of both the central atom atomic core and 
the ligand core. In formulas: 



^Hyb ^ I dr[ + ] p{r) 



\r-R,\ \r-R,\ 



with obvious meaning of the symbols. Therefore, the restor- 
ing force constant is given by 



where L,; is the appropriate component of the angular mo- 
ment. 
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Fig. 2. Atom witli axial symmetry. Only tlie charge external 
to the sphere of radius contributes to both the restoring 
force constant and the moment of inertia. 

4 Conclusion 

In a previous paper we suggested, in the framework of 
a semiclassical model, that axially symmetric atoms in 
crystals with ionic bond should exhibit a collective exci- 
tation which should be interpreted as a precession around 
the anisotropy axis of their cells. Such excitation has a 
magnetic character and a clear cut signature due to the 
dichroism of the photoabsorption cross-section. 

In the present work we studied how the above picture 
must be modified in crystals with cubic symmetry. This 
problem is interesting per se and because of its connection 
with compounds whose magnetic anisotropy is at the ori- 
gin of many technological applications. Among such com- 
pounds some might have a Spin-orbit Locking structure 
which, as already anticipated, might be investigated by 
excitations of the Scissors Modes [M]. 

The main result is that a collective mode of magnetic 
character can be expected also in such crystals, with a 
photoabsorption cross-section of the same order of mag- 
nitude as that for crystals with axially symmetric atoms. 
This cross-section, however, does not depend on the an- 
gles between the direction of the photon and the axes of 
the cell. 

The evaluation of the parameters of the model is much 
more difficult in the presence of a covalent bond, in which 
case the very existence of the Scissors Mode becomes un- 
certain. Indeed its basic assumption is that the moving 
atoms should not change their shape. If the shape changes, 
intrinsic motions of electrons will couple to the collective 
motion, with the possibility of a significant fragmenta- 
tion. In fact while in Bose-Einstein condensates the Scis- 
sors Mode fully shows its collective character [S] , in atomic 
nuclei it is split into two or more close levels. If fragmenta- 
tion is too high, the magnetic strength is distributed on a 
quasi continuum of levels and the collectivity of the state 
is altogether disrupted. The question therefore is to which 



extent semiclassical states of the type considered can exist 
in a fully quantum framework as a covalent bond. This and 
the problem of a microscopic evaluation of the moment of 
inertia are issues of conceptual interest. 
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